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ABSTRACT 

Context. Despite their profound eff'ect on the universe, the formation of massive stars and stellar clusters remains elusive. Recent 
advances in observing facilities and computing power have brought us closer to understanding this formation process. In the past 
decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, 
the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against background mid-IR 
emission, whereas dust continuum observations allow cold, dense pre- stellar-clusters to be identified anywhere. 
Aims. We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that 
are not mid-IR-dark, and to roughly characterize the level of star formation activity in dust continuum sources. 
Methods. We use Hi-GAL 70 to 500 yt/m data to identify dust continuum sources in the ^=30° and ^=59° Hi-GAL Science 
Demonstration Phase (SDP) fields, to characterize and subtract the Galactic cirrus emission, and perform pixel-by-pixel modified 
blackbody fits on cirrus- subtracted Hi-GAL sources. We utilize archival Spitzer data to indicate the level of star-forming activity in 
each pixel, from mid-IR-dark to mid-IR-bright. 

Results. We present temperature and column density maps in the Hi-GAL ^=30° and ^=59° SDP fields, as well as a robust algorithm 
for cirrus subtraction and source identification using Hi-GAL data. We report on the fraction of Hi-GAL source pixels which are 
mid-IR-dark, mid-IR-neutral, or mid-IR-bright in both fields. We find significant trends in column density and temperature between 
mid-IR-dark and mid-IR-bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on 
average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to 
star-forming, and that warmer sources are associated with more star formation tracers. Additionally, there is a trend of increasing tem- 
perature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 jum bright 
sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far- side of the Galaxy. These are cold (~ 20 K), 
high column density (N(H2) > lO^^cm"^) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show 
little to no absorption at 8 yum. These are the first inner Galaxy far-side candidate IRDCs of which the authors are aware. 

Key words. ISM: clouds - ISM: dust, extinction - ISM: evolution - stars: formation - stars: pre-main- sequence 

1, Introduction proto-clusters is complicated: they are rare, and therefore fur- 
^ ther away on average than isolated low mass star-forming re- 
5^ Massive stars play a dominant role in shaping the Universe gi^^s, proto-clusters have high column densities, meaning that 
through their immense lomzing radiation, wmds, and spectac- ^^e proto-stars are highly embedded, and massive stars evolve 
ular explosive death, yet their formation mechamsm remams j^^pidly and quickly heat and ionize their surroundings, disrupt- 
poorly understood. The dominant mode of star formation, and jj^g ^j^^jj. x^^i^i molecular cloud 
perhaps t he only mode for massive star formati on seems to be 

clustered (|Lada & Lada|2003|; Ide Wit et al.ll2005h. The definition Searches for proto-clusters or massive proto-stars are usually 
of clustered can be called into q uestion (e.g.|Bressert et al.||2010|; targeted at long- wavelengths, as dust continuum emission from 
iGieles & Portegies Zwaril2011). but the search for young, mas- these cold, dense sources peaks in the far-IR to sub-mm. Surveys 
sive star forming regions is still directed toward 'proto-clusters:' of molecular lines are another approach for identifying poten- 
cold, dense, massive molecular clumps. Direct observation of tial proto-clusters, although these surveys may be time-intensive 
and the strength of a molecular line is complicated by exci- 



* Herschel in an ESA space observatory with science instruments tation conditions and o ptical depth The d iscovery of Infrared 
provided by European-led Principal Investigator consortia and with im- Dark Clouds (IR DCs lEgan et al.1 1 19981 : IPerault et al.1 1 19961: 
portant participation by NASA. lOmont et"aDl2003h . opened a new window to viewing these cold, 
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dense potential proto-clusters in silhouette against the bright 
Galactic mid-IR background. In the past decade, compelling ev- 
idence has emerged that suggests that some IRDCs may be pre- 
cursors to niassive s tars and clusters (e.g . Rathborne et al. 2006; 



Ragan etal.] 120061: iBeuther & Sridharan .2007: Parsons et al.. 
20091: iBattersbv et al.1 l2010h Whi l e there exist many smal l 
IRDCs (e.g. iPeretto & Fullell2010l: iKauffmann & Pillaill2010h . 
the most massive on es (M ~ 10^""^ Mp, nn > 10^ cm"^, T < 
25 K jRathborne et al .' 2006: E gan et al.lll99llCarev et al.ll 19981) 
are consistent with expectations for a proto-cluster. However, the 
identification of an IRDC requires that it be on the near- side of 
a bright mid-IR background. This limits our potential to under- 
stand the Galactic distribution of potential proto-clusters. 

Far-IR and submm dust continuum surveys are a powerful 
way to identify proto-clusters throughout the Galaxy, as the cold, 
dense dust is optically thin at these wavelengths. Surveys such 
as th e Bolocam Galactic Plane Survey (BGPS, lAguirre et al.l 
at 1.1 mm, the AP EX Telescope Large Area Survey of 
the Galaxy (ATLASGAL. ISchuUer et'aDl2009l) at 870 fim, and 
now Hi-GA L (The Herschel Infrared Galactic Plane Survey, 
iMolinari et al. 2010a) from 70 to 500 yum are promising tools 
for understanding star cluster formation on a Galactic scale. 
However, the sources identified in these surveys may span a 
large range of evolutionary states, from pre- star-forming to star- 
forming, and further analysis or intercomparison may be nec- 
essary to identify the pre- star-forming regions. Since dust tem- 
peratures and column densities can be derived from the multi- 
wavelength Hi-GAL data, this data set allows for the distinction 
to be made between pre- and star-forming regions. 

For this stu dy, we utilize data from the Hi-GAL survey 
dMoliimi et al.l l2010a) from 70 to 500 yum in the ^=30° and 
^=59° SDP fields to characterize dust continuum sources. We 
investigate diff'erences in the physical properties of mid-IR-dark 
and mid-IR-bright clouds identified in the dust continuum, and 
also the physical properties of sources associated with vari- 
ous star formation tracers. We use Extended Green Objects 
(EGOs, also known as "green fuzzies' MCyganowski et al ]|200a 
[Chambers et al.l l2009 l) to trace outflows from young stars, 
CH3OH masers to trace sites of massive star formation, and 8 
and 24 i im emission to indicat e an accreting proto-star or UCHII 
Region (IBattersbv et al.ll2010h . 

This paper is organized as follows. In Section [2] we intro- 
duce the Hi-GAL observing strategy and discuss archival data 
used in our analysis. In Section [3] we describe the Galactic cirrus 
emission removal and source identification methods, the mod- 
ified blackbody fitting procedure, and how the star formation 
tracers were incorporated. Section |4] describes our results, in- 
cluding a discussion of uncertainties, the properties of the cirrus 
cloud emission, the temperature and column density maps, the 
association of Hi-GAL sources with mid-IR-dark and mid-IR- 
bright sources, and star formation tracers in Hi-GAL sources. In 
Section[5]we present five candidate IRDC-like clouds on the far- 
side of the Galaxy, and finally, in Section [6] we summarize our 
conclusions. 



2. Observations and Archival Data 

2.1. Hi-GAL 

The Herschel Infrar ed Galactic Plane Survey, Hi-GAL 
(IMolinari et al.l bOlOa '). is an Op en Time Key Project of the 
Herschel Space Observatory (Pilb ratt et al.| [2010). Hi-GAL will 
perform a 5 -band photometric survey of the Galactic Plane in a 
\b\ <r -wide strip from -70° < 1 < 70° at 70, 160, 250, 350, 



and 500 //m using the PACS (iPoglitsch et al.ll2010l) and SPIRE 
(iGriffin et al.ll2010l) imaging cameras in parallel mode. Two 2° x 
2° regions of the Hi-GAL survey w ere complet e d durin g the 
Science Demonstration Phase (SDP; Molinari et al.ll2010ab . cen- 
tered at approximately [^,b] = [30°, 0°] and [59°, 0°]. 

Data reduction was carried out using the Herschel Interactive 
Processing Environment (HIPE, Ott 2010) with custom reduc- 
tion scripts that deviated conside rably from the standard pro- 
cessms for PACS (Poditsc h et al.ll2010h . and to a lesser extent 
for SPIRE CGriffin et aL.201Qh . A detailed descr iption of the en- 
tire data reduction procedure can be found in iTraficante et al.l 
(2010). 

The zero-level off'sets in the Herschel maps were estab- 
lished by comparison with the IRAS and Planck data at com- 
parable wavelen gths, fo llowing the same procedure as de- 
scribed in Bernar det al.l (^010). We compared the Herschel- 
SPIRE and PACS data with the predictions of a model pro- 
vided by the Planck collaboration (Planck Core-Team, private 
communication) and constrained on the Planck and IRIS data 
(iMiville-Deschenes & Lag ach^ 12005 ). The model uses the all- 
sky dust temperature maps derived from the IRAS 100 jim and 
the two highest Planck frequencies to infer the average radi- 
ation field intensity for each pixel at the common resolution 
of the Planck and IRA S resolution of 5'. The Dustem model 
(Compiegne e t al.l 12010 ^ with the above value for the radiation 
field intensity was then used to predict the expected brightness 
in the Herschel- SPIRE and PACS bands, using the nearest avail- 
able Planck or IRAS band for normalization and taking into 
account the appropriate color correction in the Herschel filters. 
The predicted brightness was correlated with the observed maps 
smoothed to the 5' resolution over the region observed with 
Herschel and the offsets were derived from the zero intercept 
of the correlation. We estimate the accuracy of the off'set deter- 
mination to better than 5%. 

2.2. Archival Data 

We utilize the wealth of archival data in the Galactic Plane for 
our analysis. We use the mid-IR data taken as part of the Galactic 
Legacy Infrared Mid-Plane Su rvey Extraordinaire (G LIMPSE; 
3.6, 4.5, 5.8, and 8.0 fim, Benjamin et alj l2003h and the 
MIPSGAL survey (24 //m; Carev et al. 2009). also make use 
of the IRDC catalog of Peretto & Fuller ( 2009), the catalog of 
Extended Green Objects (Cyganowski et al. 2008, EGOs), and 
the 1.1 mm data and catalog from th e Bolocam Galacti c Plan e 
Survey (BGPS; lAguirre et al] l2010l: iRosolowskv et al.1 l2010l) . 
The Multi-Array Galactic Plane Im aging Survey (MAGPIS; 
Helfand et al. 2006; White et al. 2005), which provides compre- 
hensive radio continuum maps of the first Galactic quadrant at 
high resolution and sensitivity, is also used in our analysis. 

3. Methods 

3.1. Source Definition and Removal of the Galactic Cirrus 

The Hi-GAL data reveal a wealth of structure in the Galactic 
Plane, from cirrus clouds (Martin etal. 2010) to filaments and 
clumps, as discussed in Molinar i et al.l (12010b) . Figure[T] demon- 
strates the complicated association of mid-IR and dust contin- 
uum sources toward the Galactic Ring. In this paper, we explore 
the physical properties of the densest components of the Galactic 
Plane; the potential precursors to massive stars and clusters. In 
order to properly characterize these dense objects, a careful re- 
moval of the Galactic cirrus is required. We have explored a va- 
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Fig. 1. A three-color GLIMPSE image with logarithmic BGPS 1.1 mm dust continuum contours on our test region in the Hi-GAL 
SDP ^=30° field. Red: 8 yum , Green: 4.5 yum , and Blue: 3.6 yum. This figure demonstrates the distinction between the IRDC and dust 
continuum population. Some dust continuum sources are mid-IR-bright, some are mid-IR-dark, and some have no IR correlation. 
The object in the center is particularly well- suited as a test case for this study as it contains both mid-IR-bright and mid-IR-dark 
dust continuum sources radiating as filaments on either side of a young HII region complex. 



riety of methods for the removal of the cirrus emission and we 
briefly discuss the pros, cons, and systematic efl'ects of the dif- 
ferent subtraction methods. 

Our first step is to project all the data onto a common grid, 
with a common resolution and common units for comparison. 
We crop the images to the useable science field, convert the data 
to units of MJy/sr, Gaussian convolve to a common resolution 
(36 and 25 for the lower and higher resolution modified 
blackbody fits), and regrid to a common grid with a reasonable 
pixel size (roughly 1/4 the &fwhm)- These images are used in 
the remainder of the analysis. We include images throughout the 
paper of a "test field," (approximately, [^,b] = [29.95, -0.35] to 
[30.50, 0]) chosen for it's diversity of mid-IR-bright and mid- 
IR-dark sources. The analysis was orginally run and optimized 
on this region and then expanded to the full Hi-GAL fields. 

The first method for determining the dense clump source 
masks was to use Bolocat clump mas ks derived from the 
Bolocam Galactic Plane Survey (BGPS lAguirre et al] l2010l: 
lRosolowskv etal.1 12010^ at 1.1 mm. While these masks are a 
robust tracer of the cold, dense gas, they only extend to \b\ < 
0.5° in the majority of our science field and the sensitivity in the 
^=59° field is too poor to trace the clumps seen in Hi-GAL. 

The following methods all use the SPIRE 500 jim data 
to determine the source masks, as the cirr us decreases to- 
wards longer wavelengths (iGautier et al.ll 19921) and lPeretto et al.l 
(l2QTQh demonstrate that the SPIRE 500 //m data is well- suited 
for cirrus/dense source distinction. The second method we tried 
was based on that of Peretto et al. ( 2010), who found that toward 
IRDCs, the SPIRE 500 jim flux distribution showed a clear peak 
at low fluxes with a long tail toward higher fluxes. While this 



metho d works well when applied to a single IRDC, as it was 
used in lPeretto et al.l (l2010l) . it is not robust enough to be used to 
create a full Galactic cirrus emission map from high (\b\ > 1 ° ) to 
low (\b\ = 0° ) Galactic latitudes. We then tried a third method to 
create the dense source masks using a simple contrast map (con- 
trast = data - background) of the SPIRE 500 jim data and applied 
a cutofl'. This method improved over the previous, however, it 
sufl'ered from large negative bowls, creating source masks that 
were much smaller than the physical source sizes. In the fourth 
method, we define the background first by a second-order poly- 
nomial plane fit (along Galactic latitude) to the smoothed 500 
Jim image. This plane fit to the Galactic cirrus was then sub- 
tracted from the original data and a cutoff' was applied (deter- 
mined by eye) to define the dense source masks. This was the 
first method to do a reasonable job of identifying sources in a 
range of Galactic latitudes. We found, however, that a polyno- 
mial plane was not a very good approximation to the shape of the 
Galactic cirrus emission, and that the fit was particularly poor at 
high Galactic latitudes (\b\ ~ 0.8 - 1 ° ). 

Figure [2] depicts the final method used to identify the dense 
sources and to separate those from the cirrus cloud emission. The 
first panel (a) shows the original SPIRE 500 yum image. The first 
step (shown in panel b) is to convolve the 500 yum image with 
a Gaussian that is large enough to smooth over the sources but 
small enough to capture variations in the cirrus emission. We 
decided through trial and error that a Gaussian with a FWHM 
of 12' was a reasonable compromise. We then fit a Gaussian in 
latitude to each Galactic longitude, as shown in Figure [3] and in 
panel (c) of Figure [21 A Gaussian is a reasonable approximation 
for the variation in the Galactic cirrus across the Galactic Plane 
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and worked to identify sources at high and low Galactic latitudes 
equally well. We then subtract the Gaussian fit approximation of 
the Galactic cirrus from the original 500 jim data to achieve a 
"difference image." This "diff'erence image" is the first guess at 
the cirrus- subtracted source map. 

We apply a cutoff' to this first guess cirrus- subtracted 
source map, such that everything above the cutoff' is considered 
'source.' We determine this cutoff' by fitting a Gaussian to the 
histogram of pixel values (as shown in Figure |4]). The negative 
flux values in the distribution of this "diff'erence image" are rep- 
resentative of the random ffuctuations in the data, so we mirror 
the negative ffux values about zero, fit a Gaussian to that distri- 
bution (representative of the noise in the map) and then apply 
a cutoff' of 4.25 cr to the "diff'erence image." We calculated this 
cutoff' in quarter cr intervals from 3 to 6 cr, a range over which 
it grows smoothly. The choice of 4.25 cr was selected because 
that cutoff' best represented the sources in both the ^=30° and 
^=59° fields when inspected by eye. While the choice of cutoff' 
is important in determining the final physical properties, there is 
nothing special about this cutoff'; the properties vary smoothly 
above and below this value. 

Once the source cutoff' has been determined, the sources are 
masked out in the original 500 yum image (as shown in solid 
white in panel (d) of FigureO and we convolve (with a Gaussian 
of FWHM 12') the data outside the source masks, the cirrus 
emission, treating data in the source masks as missing. This cre- 
ates a second guess at the cirrus emission, equivalent to panel 
(b). We then repeat the process of fitting a Gaussian in latitude 
to this image (as in panel (c)), determining a source cutoff' by fit- 
ting a Gaussian to the mirrored negative ffux distribution in the 
"diff'erence image," applying these masks (panel (d)), and con- 
volving around the masks (panel (e)) to produce the next guess 
at the cirrus emission (panel (b)). We iterate on this process until 
the source masks converge. Figure [5] shows the convergence of 
the source cutoff' in both fields. We chose iteration 16 as the fi- 
nal source mask cutoff' as it was representative of the converged 
value around which the iterations varied, though the exact choice 
does not matter much, as the cutoff' varies very little after itera- 
tion 10 in both fields. In panel (d) of Figure [21 the white masks 
are the first iteration source masks and the white contour shows 
the final (16th) iteration source masks. Once the source masks 
are determined, they are applied to each wavelength image, and 
that image is convolved (with a Gaussian FWHM 12' and the 
mask pixels ignored) to create the cirrus image at that wave- 
length. Subtracting the cirrus image at each wavelength from the 
original data produces the source images used in the modified 
blackbody fits, while the smooth cirrus images are used for the 
modified blackbody fits to the cirrus emission. 

The top panel of figure [3] shows a slice of 500 jim data 
through Galactic latitude. The final convolved background 
(dashed black line) is representative of the low-lying emis- 
sion across Galactic latitude, including it's asymmetry about 
b=0° . Notice in the bottom panel that the structure is fiat across 
Galactic latitude and allows for the identification of the source 
near b=l° . Figure IH shows that as the iterations converge, we 
move ffux from negative bowls back to sources as positive fea- 
tures and the mirrored negative distribution then more closely 
resembles a Gaussian (the black line is the final Gaussian fit to 
iteration 16). 

This method has been fine-tuned by eye to reproduce the sig- 
nificant structure picked out by a human, but is entirely auto- 
mated and was run in the exact same way for both the ^=30° and 
^=59° fields, despite their vast diff'erences. We have found that 
this method faithfully identifies sources across the range of 
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Fig. 3. A slice across the Galactic Plane at 1 = 30.51° of the 
SPIRE 500 Jim image, demonstrating the cirrus emission re- 
moval. The top panel shows the SPIRE 500 yum data, with a 
Gaussian fit to the smoothed background on the first iteration 
shown in a solid gray line. The dashed gray line is the convolved 
background on the first iteration, after the source mask has been 
applied. The black solid line is the Gaussian fit to the smoothed 
background on the final iteration and the black dashed line is 
the final convolved background after the source mask has been 
applied. Note that the final background (dashed black line) fits 
the low-lying diff'use emission nicely, including it's asymmetry 
about b = 0° . The bottom panel shows the final background 
subtracted science image cut along the same Galactic longitude, 
with the final source cutoff' drawn as a dashed black line at 65 
MJy/sr. The success of the background subtraction method is 
demonstrated by the fact that the structure is ffat across Galactic 
latitude and that the cutoff' allows for the identification of the 
significant, though faint, source near b=l° . 



Galactic latitudes covered by the data and the range of source 
confusion and activity observed. The sensitivity of this method 
is dependent upon the final mask cutoff' or confusion, rather than 
the sensitivity of Hi-GAL. For the ^=30° final mask cutoff of 65 
MJy/sr, we are sensitive above the cirrus background to a col- 
umn density of N(H2) = 2.9 x 10^^ cm'^ for 20 K dust or N(H2) 
= 9.7 X 10^0 cm-2 for 40 K dust. In the ^=59° field the final 
mask cutoff' is about 24 MJy/sr, which gives a sensitivity above 
the cirrus background to a column density of N(H2) = 1.1 x 10^^ 
cm-2 for 20 K dust or N(H2) = 3.6 x 10^^ cm'^ for 40 K dust. 

3.2. Modified Blackbody Fits 

We performed pixel-by-pixel modified blackbody fits to both the 
cirrus cloud and the dense source clump emission. The dense 
source pixel-by-pixel fits were first performed on the images 
convolved to 36" resolution. Following that, the results from the 
36" resolution tests were used as input to pixel-by-pixel fits per- 
formed on the images convolved to 25" resolution. While the 
fits to the 25" resolution images formally utilize less data points, 
they allow for a more detailed comparison. 
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Fig. 2. A depiction of the cirrus subtraction method for the first iteration. Panel (a) is the original SPIRE 500 yum image (all images 
on the same linear reverse grayscale, from -10 to 900 MJy/sr). Panel (b) is the smoothed (FWHM Gaussian convolution kernel of 
12') SPIRE 500 yum image that is used to fit the Gaussian shown in panel (c). Panel (c) is then subtracted from panel (a) to produce 
a contrast image. A 4.25 cr cutoff' is then applied to the contrast image to produce the source masks shown in panel (d). Panel 
(d) is the original SPIRE 500 jum image with the sources masked out in white. This is the first iteration source masks; the final 
source boundaries are shown as white contours. Panel (d) is then convolved with a Gaussian FWHM convolution kernel of 12' with 
the masks treated as missing data to produce panel (e). Panel (e) is considered the cirrus image for the next iteration. In the next 
iteration, a Gaussian will be fit to the cirrus image, as a function of latitude at each longitude, and steps b - e will repeat. Panel (f) 
is the diff'erence image (original SPIRE 500 yum image - convolved cirrus image) that will eventually be used for greybody fitting in 
the final iteration. 



We use the modified blackbody expression in the form 
2hy^ 



c^(ekT - 1) 



-(1-^-^0 



where 



Ty =yUH2mH/^yN(H2) 



(1) 



(2) 



where /uh2 is the m ean molecular weight fo r which we adopt a 
value of JIH2 = 2.8 (iKauffmann et al.ll2008l) . mn is the mass of 
hydrogen, N(H2) is the column density, and Ky is the dust opac- 
ity. We determine the dust opacity as a continuous function of 
frequency by fitting a power-law of the form Ky = ko(^)^ to the 
tabulated lOssenkopf & HenninS (1 1 9941) dust opacities in the rel- 
evant range of frequencies. We used the lOssenkopf & HenningI 
MRN distribution with thin ice mantles that have coag- 
ulated at 10^ cm"^ for 10^ years model for dust opacity, which 
is a reasonable guess for these cold, dense clumps. We have as- 
sumed a gas to dust ratio of 100, which yielded a kq of 4.0 at vq 
of 505 GHz withayS= 1.75. 

We first fit a modified blackbody to each pixel in the cir- 
rus image (e.g. Figure [2l panel (e)) using MPFITFUN (from the 
Markwardt IDL Libarary, Markwardt 2009). We assign a cali- 
bration uncertainty of 20% to the SPIRE and PACS data points 
and use the covariance matrix returned by MPFITFUN to esti- 
mate errors in the parameters. In the modified blackbody fit, we 
leave the column density, JS, and the temperature as free parame- 
ters. The temperature was restricted to range between and 100 
K, the column density was free to range between and 100 x 



10^^ cm-^ and ft from 1 to 3. We do not use the BGPS 1.1 mm 
point in the fits because uncertainties in both the absolute cal- 
ibration and the spatial filter function prevent the BGPS point 
from significantly helping to constrain the fit at this time. 

We then fit a modified blackbody at each position within the 
source masks of the cirrus- subtracted science images (e.g. Figure 
[21 panel (f)) to the dense clumps. For these fits of dense clumps, 
we ignore the 70 yum point since the optically thin assumption 
may not be valid (rvo^^m = 1 at N(H2) = 1.2 x 10^^ cm"^ us- 
ing Equation [21). In fact, IRDCs often appear in absorption at 70 
yum indicating that in these very high column density regimes, 
the 70 yum data is not optically thin and that it can no longer 
be modelled by a sing le dust temperature. Additionally, it has 
been found (e.g. D esert et al.l lT990: ?) that a large fraction of 
the ffux at 70 yum is likely contaminated by emission from very 
small grains whose temperature ffuctuations with time are not 
representative of the equilibrium temperature of the large grains. 
Without the 70 yum point we have only four points (usually just 
on the Ra yleigh- Jeans slope) to constr ain the fit, so we fix j3 to 
1 .75 rThe lQssenkopf& HenningI 1994 fit value for an MRN dis- 
tribution with thin ice mantles that have coagulated at 10^ cm"^ 
for 10^ years; a reasonable guess for these environments), and 
leave only the temperature and the column density as free pa- 
rameters. We discuss how a value of of 1 .5 or 2 would alter the 
results in Section lTTl We acknowledge that perceived changes in 
temperature may actually be changes in JS and discuss this more 
also in Section |4J1 We are, to an extent, looking for changes in 
the dust properties (J3 or temperature) across diff'erent environ- 
ments, which we can still achieve. 
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Difference Image Flux [MJy/Sr] 

Fig. 4. The flux distribution in the difl'erence images (original 
500 yum image - convolved cirrus image, e.g. Figure [2] panel (f)) 
for 8 iterations from 1 to 20 in the ^=30° field. In each itera- 
tion, the negative flux distribution was mirrored about and a 
Gaussian was fit to the distribution to determine the character- 
istic fluctuations so that a source identification cutoff' (of 4.25 
cr) could be applied. Plotted here are the full distributions, not 
the mirrored distributions used to fit the Gaussian. As we iterate, 
points in the negative end are transferred to the positive end of 
the distribution as ffux is restored in the negative bowls around 
bright sources. The Gaussian fit to the distribution for the final 
iteration (16) is shown as the black curve. While the final dis- 
tribution is not perfectly Gaussian, the iterations show great im- 
provement, and a Gaussian that is characteristic of the random 
fluctuations in the distribution. We note that the distribution does 
not peak at zero because the average must be zero and there is a 
large positive wing in this distribution. 
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Fig. 5. The 4.25 cr mask cutofl' over 20 iterations. The cutofl' is 
from the Gaussian fit (see Figure ^ to the flux distribution of 
points in the difl'erence image (e.g Figure [2] panel (f)). Both the 
^=30° and ^=59° fields converge after several iterations, but to 
very diff'erent values. We chose iteration 16 as the final cutoff' as 
it is high enough to be assured of convergence and is representa- 
tive of the value to which the iterations converge. 

These pixel-by-pixel modified blackbody fits to the cirrus- 
subtracted dense clumps were first performed on the images con- 
volved to 36" resolution. Following that, the results from the 
36" resolution tests were used as input to pixel-by-pixel fits per- 
formed on the images convolved to 25" resolution. While the 
fits to the 25" resolution images formally utilize less data points 



(the SPIRE 500 fim point must be excluded because of resolu- 
tion), they allow for a more detailed comparison. The physical 
properties returned by the fits at diff'erent resolutions agree well. 
The final temperature and column density maps (as shown in 
Figure [6]) show the derived temperatures and column densities 
from the 25" map within the source masks, and the smoothed 
cirrus-background temperatures and column densities outside 
the source masks. 

We note here that star-forming regions contain structure on 
many scales, and that fitting a single temperature and column 
density over a region of 25" (about 0.5 pc at a distance of 4 
kpc) is a vast oversimplification. There certainly exist large vari- 
ations in the temperature and column density in these sources 
on smaller scales. In this paper, we present the beam-diluted 
average source properties on 25" scales. The column densities 
presented here are therefore characteristic of the larger- scale 
clump structure and under-estimates of the peak column den- 
sities, while the temperatures will be underestimates in very hot 
regions (hot cores) and overestimates in very cold regions (star- 
less IRDC cores). 

3.3. Star Formation Tracer Label Maps 

In order to robustly compare the physical properties determined 
in this paper with various star formation tracers, we have created 
star formation tracer label maps. These maps are on the same 
grid as the science images (temperature, column density, etc.) 
and have a binary denotion in each pixel: 1 if the star forma- 
tion tracer is present and if it is not. We create star formation 
tracer label maps for the MIPS 24 /i m emission, 8 //m emis- 
sion. Ext ended Green Objects (EGOs: [Cyganowski et al.l [20081: 
Chamber s et al.|[2b 09), and 6.7 GHz methanol maser emission 
( Pestalo zzi et al.l200 5). The MIPS 24 /im images required stitch- 
ing together using MONTAGE before the label map could be 
created. Since the ^=30° field is more populated, confused, and 
therefore has a wider range of star formation activity, we only 
perform the analysis in this field. The results of this compari- 
son can be found in Section |431 The 8 jim emission label maps, 
however, were created and analyzed for both the ^=30° and 
^=59° fields so that we could understand the relationship be- 
tween mid-IR and Hi-GAL sources in both fields (see Section 
1441) . 

The GLIMPSE 8 yum band shows emission from warm 
dust and Poly cyclic Aromatic Hydrocarbons (PAHs). Regions 
with bright 8 jim emission may contain warm, diff'use dust or 
UV- excited PAH emis sion from an Ultra-Compact HII Region 
(Batt ersby et al.|[20TQl) . The absorption of 8 jim emission by an 
Infrared Dark Cloud (IRDC) is indicative of a high column of 
cold dust obscuring the bright mid-IR background. We note that 
mid-IR-bright regions are generally indicative of active star for- 
mation, and that the peak of mid-IR lig ht can be o ff^set from the 
youngest regions of star formation (e.g. lBeuther et al. 2007). We 
denote regions of bright 8 fim emission as mid-IR-bright (mIRb), 
dark regions of absorption at 8 jim as mid-IR-dark (mIRd), and 
regions without stark emission or absorption at 8 fim as mid-IR- 
neutral (mIRn). The 8 yum emission label map has both a neg- 
ative (IRDC, mIRd), positive (mIRb), and neutral (mIRn) com- 
ponent. This map is created by first making a contrast image. 
The GLIMPSE 8 yum image is convolved with a Gaussian of 
FWHM 5' to represent the smooth, slowly vary ing background. 
5' was determined (see lPeretto & Fulle3l2009[ for more details 
on this smoothing kernel choice) to be large enough to capture 
the largest IRDCs and mid-IR-bright objects, while still small 
enough to capture the background variations. This image is then 
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Fig. 6. Temperature (left) and Column Density (right) maps in the ^=30° field test region. In these maps, the background cirrus 
emission temperature and column densities are plotted outside of the source masks, while the background- subtracted temperature 
and column densities are plotted inside the source masks. These maps are at 25" resolution and assume = L75 inside the source 
masks. 
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Fig. 7. GLIMPSE 8 yum image overlaid with contours of the clas- 
sification as mid-IR-bright (mIRb, dark blue) or mid-IR-dark 
(mIRd, cyan). The classification is only applied where a far-IR 
source exists, as denoted by the black contours. 



subtracted from the original image, which is then divided by the 
smoothed background. The resulting 'contrast map' enhances 
the contrast of mIRb and mIRd clouds above the bright mid-IR 
background. To obtain a fair comparison with the lower resolu- 
tion (~ 25") science images, this contrast image is convolved to 
the same resolution. We then apply a cutoff' such that any pixel 
with a positive contrast greater than or equal to 10% is consid- 
ered mIRb, and any pixel with a negative contrast greater than 
or equal to 5% is considered mIRd, while all others are then 



mIRn. We rely on our eyes to pick out regions of emission and 
absorption, and have selected the cutoff's that best represent a vi- 
sual interpretation of that which is dark and bright, as shown in 
Figure [71 

A portion of this label map is shown in Figure [71 This 
method does a very good job of picking out the bright and dark 
regions, and agrees reasonably well (considering the diff'erent 
m ethods, resolution used, and cutoff's) with the IRDC catalog 
of iPeretto & Fullerl (l20Q9l) . Two important biases to note with 
this method are: 1) the contrast image technique creates 'nega- 
tive bowls' around bright objects, and therefore the sizes of the 
mIRd and mIRb regions are underestimated, and 2) convolving 
the contrast image to 25" resolution causes us to miss some small 
IRDCs. Additionally, any decrements in the background will be 
denoted as being mIRd. This issue is partly resolved when we 
multiply this label map with our source masks, so that our clas- 
sifications only apply where a Hi-GAL source also exists. 

Emission at 24 yum observed using the MIPS instrument 
traces warm dust. MIPS 24 jim emission can trace warm, dif- 
fuse dust in the ISM or warm dust associated with star formation 
such as material accreting onto forming stars, the gravitational 
contraction of a young stellar object, or warm dust surrounding 
a newly formed star. The MIPS 24 yum emission label map was 
created using the same method as the 8 yum label maps, with 
the exception that we only include a positive contrast cutoff'. 
While some IRDCs remain dark at 24 yum, the negative contrast 
is very minimal when convolved to 25" resolution. Additionally, 
the highest negative contrast in the contrast image is from the 
negative bowls around bright sources, so identifying IRDCs in 
the 24 yum images using this method is not adequate. The posi- 
tive contrast threshhold for the MIPS 24 jim contrast maps was 
25%, also chosen by eye. A higher cutoff' than at 8 yum was re- 
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Table 1. EGOs Used in This Paper 



Fig. 8. GLIMPSE 3-color images (Red: 8 yum, Green: 4.5 yum, 
Blue: 3.6 jim) of EGOs in the ^=30° field. The EGOs shown in 
the right column of the 3rd and 4th row were identified as "pos- 
sible" EGOs by ICy ganowski et al.l (I2QQ 8) and the others were 
identified by two independent viewers as EGOs as described in 
Section [33l The positions are listed in Table [TJ 



quired to select the 24 jim sources above the bright background. 
These cutoff's greatly aff'ect the resulting trends. See Section |4J1 
for a discussion of the eff'ects of cutoff's on source identification. 

We include the presence of E xtended Green Objects 
(EGOs, also ca lled "green fuzzies" 'Cyganow ski et al.l I2QQ8I: 
[Chambers et al.l '2009) as a star formation tracer. EGOs are 
regions of enhanced and extended 4.5 yum emission and are 
thoug h t to be indicative of shocks in outffows (ICyganowski et al.l 
120091) . ICvganowski et al.l (l2008b catalogued EGOs throughout 
the Galaxy in the GLIMPSE fields. There are two "possible" 
EGOs in the ^=30° field whi ch we includ e in our EGO la- 
bel map. We supplemented the Cyganowski et al.l (120081) catalog 
with an independent visual search of the ^=30° field for EGOs 
that may have been missed. An independent, cross-checked 
search identified eight further candidates, which we include in 
our analysis. A list of the locations of the EGOs used in this 
analysis is given in Tabled and they are shown in Figure [8l For 
each EGO, a box was drawn around the region that contains the 



EGO Name 


R.A. (J2000) 


Decl. (J2000) 


1 X b Size(") 


Ref.^ 


G29.23-0.05 


18:44:51.4 


-03:18:41.7 


14x 12 


1 


G29.25-0.69 


18:47:10.7 


-03:35:42.0 


14x 14 


1 


G29.84-0.47 


18:47:28.8 


-02:58:03.0 


9x9 


2 


G29.91-a81 


18:48:47.6 


-03:03:31.1 


8x8 


1 


G29.96-0.79 


18:48:49.8 


-03:00:28.7 


32x32 


2 


G30.40-0.30 


18:47:52.6 


-02:23:12.2 


17x 17 


1 


G30.42-0.23 


18:47:40.8 


-02:20:31.6 


12x 12 


1 


G30.60+0.18 


18:46:33.7 


-01:59:30.3 


14x 14 


1 


030.79+0.20 


18:46:48.2 


-01:48:54.4 


12x 12 


1 


G30.90-0.01 


18:47:45.9 


-01:49:00.0 


14x7 


1 



Notes. 1 - This w ork (see Section [33] and Figure [U, 2 - 
ICvganowski et all ( 120081) 



EGO, and each pixel in this box was given a ffag of 1, while for 
pixels outside the box the ffag is 0. 

Finally, we employ the presence or absence of 6.7 GHz Class 
II CH3OH masers as a s tar formation tracer. W e use the com- 
plete catalog compiled by'Pestalozzi et al.l (12005*) to identify un- 
biased searches for 6.7 GHz methanol masers in the ^=30° ffeld. 
While there have been numerous targeted searches for CH3OH 
masers, we utilize th e unbias ed Galactic Plane searches of 
ISzvmczak et"aD (l2002h and Elli ngsen et al.l (Il996h . on the 32m 
Torun Radio Telescope (FWHM of 5.5') and the University of 
Tasmania 26m Radio Telescope (FWHM of 7'), respectively. 
6.7 GHz methanol masers are foun d to be exclusive ly associ- 
ated with massive star formation CMinier et al.ll2003l) . and are 
often off'set from the r adio continuum emission indicative of 
an UCHII Region (e.g. IWalsh et"an Il998l: iMinier et al.l l2QQlh 
and are thought to represent an earlier stage of massive star 
formation. We chose to in clude only the unbiase d catalogs of 
lEllingsen et"aD 11996') and ISzvmczak et al.l (l2002h so as not to 
choose mIRb regions preferentially. These catalogs together give 
22 methanol maser regions, which we use to create a label 
map where the "size" of the methanol maser regions is given 
by 2x the positional accuracy R MS of the observations (30" 
and 36" for'Szvm czak et al.ll2002HEllingsen et al.lll996[ respec- 
tively). Methanol maser emission comes from very small ar- 
eas on the sky (e.g. lWalsh et"aDll998l:lMi"nier et al .' 2001). how- 
ever, methanol m asers are often clustered (ISzymczak et al. 2002i: 
lEllingsen|[T996l find multiple maser spots towards the majority 
of sources), so these "source label sizes" are meant to be repre- 
sentative of the extent of the star- forming region. 



4. Results 

We discuss here the results of the modified blackbody fits, the 
properties of the Galactic cirrus emission, a comparison of mIRb 
and mIRd pixels, and the association with star formation trac- 
ers. We discuss what defines the populations of mIRb and mIRd 
sources, and the objects which have no mid-IR association. 
Before the discussion of the results, we first mention the many 
caveats and uncertainties that accompany these findings. 

4.1. A Remark on Caveats and Uncertainties 

A major uncertainty in this analysis is the modified black- 
body fits. Firstly, star-forming regions show structure on many 
scales, and so fitting a single column density and temperature 
to any point does not adequately represent the whole region; 
rather it serves as a way to categorize and diff'erentiate between 
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Fig. 9. Top Left: GLIMPSE 8 //m image in our ^=30° test region. Right: Modified blackbody fits to the points marked on the filament 
in the top left panel. Since the 70 jim point is not included in these fits (see Section [T2l) it is plotted as an X while the other points 
are shown with their 20% calibration error bars. The errors quoted in temperature and column density are the formal fit errors. 
As discussed in Section \4~T\ these are not necessarily representative of the true errors. Bottom Left: The temperature and column 
densities derived along the filament, with 10 points marked. Note the inverse relation between temperature and column density in 
this example, especially from points 8 to 10. 



bulk, large-scale physical properties of regions. Additionally, 
this means that the quantities derived are likely to be highly 
beam-diluted. The true peak column density in any given region 
is almost certainly higher than the beam-averaged column den- 
sity reported here. Likewise, the beam-averaged temperatures re- 
ported for cold objects are overestimates of the true temperature 
minimum, while the beam-averaged temperatures reported for 
hot objects is an under-estimate of the peak. These caveats ex- 
plain the somewhat small dynamic range observed in both tem- 
perature and column density. 

In addition, the parameters derived from the modified black- 
body fits themselves are rather uncertain. For most of the 
sources, the temperature is higher than 18 K, so all the points 
used in the fit lie near the peak of the SED or on the Rayleigh- 
Jeans slope. This equates to large uncertainties in the fits, es- 
pecially at high temperatures. The median errors in the derived 
temperature are about 4 and 3 K (for the ^=30° and ^=59° fields, 
respectively), while the median errors in the derived column 
density are about 0.12 and 0.04 x 10^^ cm'^ (in ^=30° and 
^=59° fields, respectively). The error is a strong function of 
the parameter value, with a median temperature error of 13% 
and 18% and median column density error of 44% and 48% 
(in the ^=59° and ^=30° field, respectively). These uncertain- 
ties are based on the 20% calibration uncertainty assigned to the 
Hi-GAL fluxes. Additionally, the derived quantities are highly 



dependent not only on the background that is subtracted, but 
also on possible departures between the gain scale of the PACS 
and the SPIRE data. Improvements to the general issue of the 
Planck-Herschel cross calibration with respect to the accuracy 
of the ofl'set determination as well as to the gain values will be 
adressed in the future. 

The assumption of a fixed JS can be considered to contribute 
the bulk of the uncertainty in the temperat ure determination. T he 
degeneracy of /3-T has been discussed in Para dis et"aD ([2010'). If 
we adopt a value of = 1.5 (instead of 1.75), the temperatures 
are about 4K higher on average and the column densities are 
about the same. If we adopt a value ofj3 = 2.0 (instead of 1.75), 
the temperatures are about 3K lower on average, and again, the 
column densities are about the same. With only four data points 
to constrain the fits, fixing /3 is nearly essential, with the choice of 
1.75 being a reasonable guess for these environments. Variations 
in JS from 1.75 to 1.5 or 2 correspond to changes in temperature 
of about 3 or 4 K. There is the possibility that systematic changes 
in 13, for example a systematic de crease in toward the center o f 
IRDCs due to dust coagulation (lOssenkopf & Henninglll994) . 
could systematically change our results. In the dark centers of 
IRDCs, however, a variation in JS only changes the temperature 
by about 2 K (2 K hotter for = 1 .5, and 2 K colder for ft = 2.0), 
so while this is an important eff'ect to consider, it is unlikely to 
negate the observed trends in temperature. 
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Fig. 10. Top: Temperature (on a linear scale from to 50 K) and Bottom: Column Density (on a log scale from N(H2)= to 5 x 10^^ 
cm"^) in the ^=30° (left) and ^=59° (right) fields. In these maps, the background cirrus emission temperature and column densities 
are plotted outside of the source masks, while the background- subtracted temperature and column densities are plotted inside the 
source masks. These maps are at 25" resolution and assume p = 1.75 inside the source masks. 



The association with star formation tracers is also bridled 
with uncertainty. Each of the tracers included are not perfectly 
represented with our methods. The 8 and 24 yum label maps are 
entirely dependent on the cutoff's chosen. While we have done 
our best to choose reasonable cutoff's, a change in these cutoff's 
would have a significant eff'ect on the properties inferred. The as- 
sociation with EGOs suff'ers from a similar bias. We have made 
a careful attempt to identify probable EGO candidates, but the 
objects themselves are somewhat ill-defined, so again, we resort 
to the by-eye classification. The association with CH3OH masers 
is robust in that it is based on uniform sensitivity blind surveys, 
however, these surveys suff'er from large beam sizes and confu- 
sion which consequently dilute their true physical properties. All 
of these tracers also suff'er from the inherent scatter that comes 
along with varying levels of extinction and diff'erent distances. 

Despite these caveats, the trends observed over many pixels 
remain robust. The pixel-by-pixel comparison is an important 
technique for achieving the best possible associations with the 
given data, and therefore, the most robust statistics. The object 



in the center of Figure[9]highfights the importance of comparing 
values on a pixel-by-pixel basis, rather than denoting the entire 
object as either mIRb or mIRd. 

4.2. Properties of the Cirrus Cloud Emission 

While we primarily consider the cirrus cloud emission as "back- 
ground" in this paper, we report here some useful physical prop- 
erties derived in this analysis. These maps were smoothed to 
12' (&fwhm) resolution, so we report only the average prop- 
erties in Table [21 These data suggest that dust in diff'use clouds 
in the molecular ring (^=30° field) is warmer, has higher col- 
umn densities, and a steeper spectral index (fi) than dust in dif- 
fuse clouds toward the mid-outer Galaxy, outside the molecu- 
lar ring (^=59° field). In the cirrus background determination, 
we fit a Gaussian of the form: g = B -\- Ae^^~^^^ across 
Galactic latitude at each Galactic longitude. We report here the 
Gaussian fit to the median value across Galactic longitudes. In 
the ^=30° field, the best fit parameters are A = 285 MJy/sr, B = 
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Fig. 11. Normalized temperature and column density histograms of source pixels in the Hi-GAL fields separated by their 8 jim asso- 
ciation; mid-IR-dark (mIRd), mid-IR-neutral (mIRn), or mid-IR-bright (mIRb). The values printed on the plots are the medians of 
the distributions (also seen in Table [3]). The mIRd pixels tend to have the lowest temperature and highest column density, while the 
mIRb pixels have the highest temperature and lowest column density, and the mIRn pixels fall in the middle. A simple K-S test (see 
Section l44l) shows that the hypothesis that any of these populations (in temperature or column density, but not with any cutoff's, so 
the top panels only) are drawn from the same distribution can be ruled out with at least 99.7% confidence. The average temperature 
and column density is lower in the ^=59° field than in the ^=30° field. Left: Normalized temperature histograms in the ^=30° (left) 
and ^=59° fields (right). The top panels are the temperature histograms of all source pixels in the 25" resolution images. The bottom 
panels are the temperature histograms of all source pixels above a column density cutoff' (N(H2) > lO^^cm"^ in the ^=30° field, 
N(H2) > 0.3 X 10^^ cm"^ in the ^=59° field). Introducing a column density cutoff' decreases the average temperature, especially for 
the mIRb pixels. Right: Normalized percentage logarithmic column density histograms in the ^=30° (left) and ^=59° fields (right). 
The top panels are the column density histograms of all source pixels in the 25" resolution images. The bottom panels are the 
column density histograms of all source pixels with a temperature less than 25 K. These low temperature pixels have higher column 
densities on average, especially the mIRb ones. 



Table 2. Smoothed Cirrus Emission Properties 



Hi-GAL Field 




N(H2)^ 


Temperature^ 






(X 10^2 cm-2) 


(K) 


^=30° 


1.7 + 0.2 


0.7 + 0.3 


23 + 1 


^=59° 


1.5 + 0.2 


0.3 + 0.1 


21 + 1 



Notes. The values given are the median of all valid values in the 
map, and the uncertainty quoted is the standard deviation. 



38 MJy/sr, bO = -0.1° , and cr = 0.4° , while in the ^=59° field, 
the best fit parameters are A = 53 MJy/sr, B = 29 MJy/sr, bO 
= 0.1 ° , and cr = 0.4° . Previous studies (e.g. Rosolowsk v et al.l 
l20TQh have shown that the mean Galactic latitude where emis- 
sion peaks depends strongly on the Galactic latitude observed, 
but that the mean in the First Galactic Quadrant is <b>^^ - 
0.1° . The ^=30° field also shows an average Galactic latitude 
of about -0.1° , while the ^=59° field is higher than the average 
atbO=0.1° . 

4.3. The Temperature and Column Density Maps 

We present full temperature and column density maps of the Hi- 
GAL ^=30° and ^=59° SDP fields. These maps are shown in 
Figure [TOl and are available to download as FITS files online. 
The source mask label maps, the source mask temperature and 
column density maps, the error maps of these quantities as well 
as the star formation tracer label maps are all available for down- 
load. The maps presented in Figure [TOl are displayed such that 
the values inside the source masks represent the background- 



subtracted fit values, while the values outside the source mask 
are the fits to the background itself. 

The uncertainties in these maps are discussed in Section |4J1 
In addition to those, we found that in the ^=30° field, our fits 
were returning unsensible parameters at high Galactic latitudes 
(\b\ > 0.8° ). This is due to the imperfect calibration, especially at 
these low levels, so the relative ffuxes were no longer physically 
meaningful. We made several attempts at correcting or properly 
ignoring those points, with little success, and therefore, recom- 
mend that the maps in the ^=30° field only be used within \b\ < 
0.8° . Additionally, pixels near the edges of the source masks 
were just barely above the background, and can produce unphys- 
ical fits, so be cautious of any pixel that is right on the edge of 
the source mask. 

The column density follows the far-IR ffux closely, while the 
temperature is quite varied, and in many cases, inversely corre- 
lated with the column density. The median temperatures for all 
the pixels in the source masks are 26 and 20 K, while the me- 
dian column densities are 0.25 and 0.10 x 10^^ cm"^ (for the 
^=30° and ^=59° field, respectively). The highest column den- 
sity points are found in W43, the large complex near ^=30.75° , 
b=-0.05° (Bally et al. 2010), where the bright millimeter points, 
MMl - MM4, have beam- averaged column densities of around 
10^^ cm-^ 

4.4. The association of Far-IR Clumps with mIRb and mIRd 
sources 

The association of Hi-GAL sources identified in the Far-IR with 
mid-IR sources is of interest, as that gives some indication of 
their star-forming activity. We find that Hi-GAL sources span 
the range of pre- to star-forming and that there exist significant 
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Fig. 12. Left: A column density versus temperature plot comparing mIRn (8 jim neutral, magenta), mIRd (8 yum dark, blue), and 
mIRb (8 yum bright, orange) pixels within the ^=30° source masks, with the median of each value plotted as a diamond on top. 
The mIRd points occupy the coldest, highest column density portion of the plot, while the mIRb points occupy a range of column 
densities. The mIRb points average a lower column density and a higher temperature region of the plot than the mIRd points. 
The mIRn points average a lower column density than mIRd points, and a slightly higher temperature, but otherwise match that 
distribution somewhat. The mIRn points at the extreme left of the plot are due to the imperfect definition of the 8 jim masks. Where 
mIRb and mIRd regions are adjacent, the pixels between are denoted as mIRn since the two adjacent regions are blurred. This causes 
some few pixels to be denoted as mIRn which may truly be mIRd or mIRb. Right: A similar plot as the left, but with all points 
(gray), points associated with a CH3OH maser (orange) or EGO (green), and 24 [im bright points (magenta). We note that EGOs 
occupy a mid-range temperature and the highest column density portion of this plot, while the CH3OH maser points occupy a wide 
mid-range of temperatures and a high range of column densities. The 24 yum points are associated with higher temperature points 
and relatively low column density. 

Table 3. Source Properties by Their Mid-IR Association 



Hi-GAL 


Mid-IR 


Temperature^ 


High Column Temperature"^ 


Column Density^ 


Low Temperature Column"^ 


Field 


Classification 


(K) 


(K) 


(X 10^2 cm-2) 


(X 10^2 cm-2) 


^=30° 


Dark (mIRd) 


21 


19 


0.50 


0.57 




Neutral (mIRn) 


25 


23 


0.25 


0.32 




Bright (mIRb) 


34 


31 


0.23 


0.49 


^=59° 


Dark (mIRd) 


16 


14 


0.22 


0.23 




Neutral (mIRn) 


20 


15 


0.11 


0.12 




Bright (mIRb) 


30 


24 


0.08 


0.15 



Notes. The values given are the median of all valid values in the map. *^2) The median temperature of all valid points in the map with N(H2) > 
10^2 cm"2 in the ^=30° field and N(H2) > 0.3 x lO^^ cm"^ in the ^=59° field (to account for the fact that column densities in the ^=59° field are 
about 1/3 of that in the ^=30° field). The median column density of all valid points in the map with T < 25 K. 



trends in both temperature and column density between these 
populations. The growing interest in IRDCs as potential precur- 
sors to stellar clusters emphasizes the importance of identifying 
IRDC-like sources throughout the Galaxy in an unbiased man- 
ner, not just where they exist in front of a bright mid-IR back- 
ground. A survey such as Hi-GAL is essential to identify IRDC- 
like objects throughout the Galaxy while a thorough analysis is 



necessary for understanding the trends in Hi-GAL sources that 
are pre- versus star-forming. 

We find that mIRd pixels are characterized by colder temper- 
atures (colder by more than 10 K) and higher column densities 
(about a factor of two higher) than mIRb pixels (see Table [3]). 
The mIRn pixels generally have a column density similar to that 
of mIRb pixels, and a temperature midway between mIRd and 
mIRn pixels (see Figure [TT]). This trend is also apparent in the 
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Fig. 14. A comparison of star formation tracers in the ^=30° field. The identification of these tracers is discussed in Section [33l Left: 
Log histogram of the temperature distributions for all source pixels identified as mIRd (8 jjm dark, black), associated with EGOs 
(green) or a CH3OH maser (orange), mIRb (8 yum bright, blue), or 24 jim bright (magenta). The median temperature is reported in 
the legend and with the dashed vertical line for each tracer. If one assumes that a source will monotonically increase in temperature 
with time, this diagram gives a possible indication of the observed evolution of sources in this field. Right: A normalized temperature 
histogram for all pixels in the ^=30° source masks containing (blue), 1 (green), 2 (orange), and 3 (magenta) of the star formation 
tracers shown in the left panel. The pixels containing 4 star formation tracers are rare, follow the distribution of the EGOs, and have 
low number statistics and therefore are not plotted. The median temperature is again reported in the legend and with dotted vertical 
lines. This plot shows that warmer pixels are associated with more star formation tracers. 



left panel of Figure [TJ] Figure [TT] shows that introducing a col- 
umn density cutoff' to the temperature distribution decreases the 
average temperature, especially for mIRb pixels. If, on the other 
hand, we introduce a temperature cutoff' (T < 25 K) to the col- 
umn density distributions, the average column density increases. 

A simple K-S test rules out the possibility that any combi- 
nation of the mIRd, mIRn, or mIRb populations in temperature 
or column density are drawn from the same distribution with 
at least 99.7% confidence. For the K-S test, we calculated the 
eff'ective number of independent points in each distribution by 
dividing the number of pixels in each distribution by the number 
of p ixels per beam. Considering the errors discussed in Section 
14.11 these trends are significant and not likely to be a by-product 
of systematic errors. Hi-GAL is sensitive to both cold, high col- 
umn density, likely pre- stellar sources and warmer, more diff'use 
star-forming regions. We can distinguish between these using the 
temperature and column densities derived simply from Hi-GAL 
or through comparison with other tracers. 

In Figure [131 we plot the fraction of mIRd, mIRn, and mIRb 
pixels above the column density cutoff' plotted on the x-axis. In 
the ^=30° field, the mIRn fraction drops quickly as the mIRd 
fraction rises, while in the ^=59° field, mIRb fraction rises 
steadily as the mIRn fraction drops. While assigning an overall 



percentage to each of the categories is somewhat subjective, be- 
cause of the cutoff'-dependent variations, we estimate that in the 
^=30° field about 55% of the Hi-GAL identified source pixels 
are mIRd, 20% are mIRn, and 25% are mIRb. In the ^=59° field, 
we estimate that about 20% of the Hi-GAL identified source 
pixels are mIRd, 40% are mIRn, and 40% are mIRb. The fact 
that the fraction of mIRd pixels in the ^=59° field is so much 
lower than that in the ^=30° field could be an artifact of the rel- 
atively sparse mid-IR background emission in the ^=59° field, 
or it could very well be that there is a lower fraction of cold, 
high-column density clouds in the ^=59° field. 

4.5. Star Formation Tracers in Hi-GAL Sources 

We compare the temperature and column density maps with 
the star formation tracer label maps (discussed in Section [33l) . 
While we evaluate the evol ution of sources by their temperature 
alone. iBeuther et al] (|201Q|) analyze the evolution of four sources 
with changes over the entire SED. We find that the more star for- 
mation tracers associated with a source, the higher the temper- 
ature (see Figure [141 right). This trend is not surprising: a star- 
forming clump will be warmer than a pre- star-forming clump. 
It is reasonable to assume that as a clump evolves, the temper- 
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Fig. 13. Fraction of pixels in each mid-IR categorization (mIRd, 
mIRn, mIRb) as a function of column density cutoff for the 
^=30° (top) and ^=59° fields. In the ^=30° field, as the column 
density cutoff' is increased, the mIRd pixels become much more 
prevalent, replacing the common, but low-significance mIRn 
pixels. In the ^=59° field, the mIRn contribution again lowers, 
while in this case, the mIRb pixels become more prevalent. In 
both plots, the sharp cutoff' near 4 and 1 indicated where the 
number of pixels included in this analysis goes below 1000. 
We conclude that in the ^=30° field, about 55% of the far-IR 
sources are mIRd, 20% are mIRn, and 25% are mIRb, while 
in the ^=59° field, about 20% are mIRd, 40% are mIRn, and 
40% are mIRb. It is evident from these plots, however, that these 
numbers are somewhat subjective and can vary depending on the 
cutoff' selected. 



ature will increase monotonically. Therefore, we compare the 
individual star formation tracers with temperature to see if there 
is any indication about which star formation tracer may turn on 
first, and what the relative lifetimes might be. There are a lot of 
uncertainties in this analysis (discussed in 14. II) and we are com- 
paring beam-averaged physical parameters, so there is quite a lot 
of scatter. 

In the left panel of Figure [141 we see a progression of star 
formation tracers with temperature, with a huge amount of over- 
lap. The median temperatures of each population (denoted with 
a dotted line, and stated in the legend) delineate a sequence 
from cold pre-stellar clumps to warm star-forming clumps, from 
mIRd, to outffow/maser sources, to sources which are bright at 
8 and 24 fim. However, the overall distributions are very wide, 
possibly indicating long lifetimes, but more likely an indication 



of large errors in both the beam-averaged physical properties and 
the assignment of star formation tracers. 

Figure[T2] shows these trends in both temperature and column 
density space. The left panel shows that mIRd pixels are gener- 
ally colder and have higher column densities than mIRb pixels. 
The right panel shows that all the star formation tracer pixels are 
warmer than the mIRd population, with CH3OH masers in the 
mid-range of temperatures, and 24 and 8 jim bright pixels oc- 
cupying the highest temperature range. EGOs are characterized 
by a similar temperature as CH3OH masers, but a slightly higher 
column density. Interestingly, the EGOs seem to occupy two dis- 
tinct regions on the plot in Figure [121 one higher and one lower 
column density, indicating that EGOs are present in a variety of 
environments. It is not surprising that the EGOs represent the 
highest column density population, as they are the most local- 
ized star formation tracer we have used. 

We perform simple K-S tests to determine the likelihood that 
any combination of the temperature or column density distribu- 
tions of star formation tracers (pixels which are mIRd, mIRb, 
24 yum bright, or contain EGOs or masers) are drawn from the 
same distribution. For the K-S test, we calculated the eff'ective 
number of independent points in each distribution by dividing 
the number of pixels in each distribution by the number of pix- 
els per beam (about 23 in the 25" resolution images used for 
all tracers except the masers, which have an accuracy of about 
33", or 41 pixels). Since EGOs are defined locally, each pixel is 
considered independent. When we compare the temperature dis- 
tributions of all combinations of the mIRd, mIRb, 24 jim bright, 
EGO and maser populations, we can rule out with at least 99.7% 
confidence that any are drawn from the same distribution, ex- 
cept the 24 and 8 jim bright points, which have about a 55% 
probability of being drawn from the same distribution, and the 
EGO and maser populations which have about a 0.5% chance of 
being drawn from the same distribution. When we compare the 
column density distributions of all combinations of the mIRd, 
mIRb, 24 Jim bright, EGO and maser populations, we can rule at 
with at least 99.7% confidence that any are drawn from the same 
distribution, except the EGO and maser populations, which have 
about a 55% chance of being drawn from the same distribution. 
The fact that EGOs and CH3OH masers may be tracing the same 
population of sources is encouraging, as they are both supposed 
to trace young, massive outffows. We would also have expected 
that pixels which are bright at 8 or 24 jim should be tracing sim- 
ilar environments, as both can be excited by the heat or UV light 
from a young, accreting protostar. 

We see interesting trends in Hi-GAL clumps between tem- 
perature and column density and star formation tracers. We find 
that 8 yum dark pixels are, on average, the coldest, followed by 
pixels containting EGOs and/or CH3OH masers, then 8 and 24 
Jim bright points. If we assume that as a clump evolves, it will 
monotonically increase in temperature with time, then this could 
suggest a possible sequence of tracers. However, the scatter is 
still too large to suggest a definitive sequence or lifetimes. We 
should expect that one of the first stages will be the formation of 
an outffow, detectable by an EGO or CH3OH maser. Following 
that, the protostar will continue to heat its surroundings and 
light up at 24 jim. A massive protostar may form an UCHII re- 
gion while still accreting, whose UV light woul d excit e PAH s 
in the 8 yum band of GLIIMPSE. Unlike Battersbv et al.l (l201Ql) . 
we find that 24 jum sources light up around the same time as 
8 yum sources. This is almost certainly due to the use of a high, 
automated cut off' to determine the e xistence or absence of a 24 
yum source. In lBattersbv etal.1 dMTot) . UCHII regions were as- 
sociated with very bright (> 1 Jy) 24 yum point sources, which 
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Fig. 15. A possible IRDC-like source on the far-side of the 
Galaxy. This source is cold (~ 20 K), has high column den- 
sity (-3x10^^ cm"^) and is amidst the bright mid-IR Galactic 
Plane at (^,b) = (30.01, -0.27). By these measures, it should be 
an IRDC, but the 8 yum image (top right) shows only a mod- 
erate decrement and the near- side extinction column density is 
almost an order of magnitude less than expected from the far-IR 
column density. This source breaks up into a Northern (^,b) = 
(30.003,-0.266)and Southern (^b) = (30.009,-0.274) component 
at the higher resolution of the PACS 70 jim band. This source is 
a great candidate for further study. 

may be the population of bright sources we are picking out 
here. We note that even though most or all UCHII regions are 
mIRb, certai nly not all mIRb sources are UCHII regions (e.g. 
iMottram etal. 2010). The fact that the population of mIRb and 
24 fim bright are so similar is likely an artifact of the cutoffs 
chosen. The two should be associated for bright sources, which 
is the trend we see in this paper, but for dimmer sources, which 
are excluded by the cutoff, we do not know the association. The 
detection of an EGO or 24 yum point source depends on the ex- 
inction and the background, and the sensitivity for detection will 
decrease with distance. As discussed, there are many caveats and 
systematic errors to consider, however, the trends are intriguing. 

5. IRDC-like Source Candidates on the Far-Side of 
the Galaxy 

We present a list of candidate IRDCs on the far-side of the 
Galaxy in Table |4j These sources are cold (~ 20 K), high col- 
umn density (> 1 x 10^^ cm"^) objects identified in Hi-GAL 
with very weak absorption at 8 jim. These are objects that, 
were they in front of the bright mid-IR background, should be 
IRDCs, but are not. To identify these sources, we calculated an 
extinction column density map at 8 yum, using the method of 
IPeretto & Fulled (12009 ) and compared this with the Hi-GAL col- 
umn density map. Any sources with a Hi-GAL column density 
> 1 X 10^^ cm"^ and an extinction column density several (3-4) 
factors lower was considered. We performed an examination by 
eye simply to check that these are plausible candidates. The ex- 
tinction must be weak, the region cannot be mIRb, and it needs 
to be cold (< 30 K). In other words, these are sources that, if on 



the near- side of the Galaxy, should be dark at 8 yum, yet they are 
not. None of these sources are identified as IRDCs by the cata- 
log of Peretto & Fuller ( 2009). Sources at high Galactic latitudes 
(\b\>0.5° ) and in the ^=59° field were not included because the 
sparse background adds additional uncertainty to the extinction 
column density estimate. 

We discuss G030.01-0.27, a candidate far-side IRDC, in 
more detail here as an example of this class of objects. G030.01- 
0.27 is a high-column density (~ 3 x 10^^ cm'^), cold (~ 20 
K) dust continuum source in the ^=30° field near the Galactic 
mid-plane that is neither 8 yum dark nor bright (see Figure [TSl). 
There is a slight decrement at 8 jim that is nearly consistent 
with background fluctu ations as this source w as not detected 
in the IRDC catalogs of iPeretto & FuUed (l2009h or lSimon et all 
(12006 ). At 70 yum, this source breaks up into two, a Southern 
source at (^,b) = (30.009,-0.274) and Northern source at (^,b) = 
(30.003,-0.266). The Southern source appears to be very small 
and faint at 8 yum, while the Northern source is a faint, slightly 
extended source at 8 yum. The Southern source has bright 24 and 
70 jjm emission. The Northern source is also seen in emission 
at these wavelengths, thou gh le ss brightly. Comparison with the 
MAGPIS (IHelfand et al.ll2006l: I White et al. 1 120051) 20 cm maps 
show no emission in the radio continuum coincident with either 
source, while the MAGPIS 6 cm map shows a very faint object 
associated with the Southern source. 

The faint decrement at 8 yum can be morphologically 
matched with a molecular cloud in ^^CO (iJackson et al.l 120061 
from the Galactic Ring Survey) at about 104 kms"^ Integrating 
from 100 to 10 8 kms~ \ assuming 20 K, and using Eq. 17 from 
iBattersbv etaP dlOlO") the column density is about 2.5 xlO^^ 
cm"^, consistent with the Hi-GAL column density estimate. 
Using the rotation curve o f iReid et al.l (l2009l) . the near- side dis- 
tance of this ^^CO cloud is 5.8 kpc and the far- side distance is 
8.7 kpc. Usin g the extinction mass determination technique of 
iButler & Taiil 12009) as applied in IBattersbv et al.l (^2010), the 8 
fim extinction column density of the cloud is 0.5 xlO^^ cm"^ as- 
suming it is on the near- side of the Galaxy or 1 xlO^^ cm"^ as- 
suming it is on the far- side of the Galaxy. It is not expected that 
the extinction column density estimation method will be very ro- 
bust for clouds on the far- side of the Galaxy, but this very slight 
extinction is at least consistent with the cloud being on the far- 
side of the Galax y. Additio nally, this cloud shows no obvious HI 
self-absorption OStil et al.l l2006). While none of this evidence is 
conclusive, it all seems to indicate that this cloud is on the far- 
side of the Galaxy. Another possibility, however, would be that 
this cloud is simply behind the majority of the diffuse mid-IR 
emission, and still on the near-side of the Galaxy. Either way, 
it is an interesting source that was otherwise missed as such in 
previous study. 

The fact that the mid-IR background is bright at the loca- 
tions of these sources, yet the 8 yum extinction column density 
is extremely low is a strong indication that these sources are be- 
hind the bright mid-IR background. In fact, the Hi-GAL distance 
analysis by ? finds that all the sources in Tableware either at the 
far or tangent distance (> 7 kpc) except for G029. 28-0.33, which 
is on the near-side at about 6 kpc, so still likely behind most of 
the bright mid-IR emission. 

With their high far-IR column density and low temperature, 
these sources are easily categorized (based on the left panel 
of Figure [T2| as IRDC-like. These are the first candidate far- 
side IRDCs of which the authors are aware. The potential outer 
Galaxy IRDC identified by Fries wijk et al. (2007) is similar in 
that they are both IRDC-like objects identified independent of 
8 yum absorption. However, these sources are not in the outer 
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Table 4. Far-Side IRDC Candidates 









Hi-GAL 


Extinction^ 






1 


b 


Column 


Column 


Temperature 


Source Name 


n 


n 


(xl022cm-2) 


(xl022cm-2) 


(K) 


G030.60+0.18 


30.602 


+0.177 


6.2 


1.1 


22 


G030.34-0.il 


30.345 


-0.114 


1.7 


0.15 


19 


GO3O.OI-O.272 


30.006 


-0.270 


3.1 


0.5 


21 


G029.31-0.05 


29.312 


-0.051 


1.0 


0.06 


13 


G029.28-0.33 


29.284 


-0.331 


1.6 


0.02 


20 



Notes. Not knowing the distances to all these sources, we use the lPeretto & FulleJ (l2009h extinction column density estimation method, which 
basically assumes the cloud is on the near-side of the Galaxy This is the source that is discussed in the text in Section[5] 



Galaxy, but on the far- side of the inner Galaxy. These candidates 
are just a few of potentially many more IRDC-like sources which 
remain to be uncovered in a study like this over the Galactic 
Plane. 

6. Conclusions 

We have performed cirrus- subtracted pixel-by-pixel modified 
blackbody fits to the Hi-GAL ^=30° and ^=59° SDP fields. The 
source identification and cirrus- subtraction routines are robust 
and can be applied to Hi-GAL data throughout the Galactic 
Plane. We present temperature and column density maps of the 
dense clumps in these fields and cirrus column density, tempera- 
ture, and 13 maps. We find that the cirrus cloud emission is char- 
acterized hy/3= 1.7, N(H2) = 0.7 x lO^^cm-^, and T = 23 K in 
the ^=30° field and /3 = 1.5, N(H2) = 0.3 x lO^^cm'^ and T = 
21 Kin the ^=59° field. 

We also characterize each pixel as mid-IR-bright (mIRb), 
mid-IR-dark (mIRd), or mid-IR-neutral (mIRn), based on the 
contrast at 8 jim. The association of Hi-GAL sources identified 
in the far-IR with mid-IR sources is of interest, as the far-IR 
sources span the range of pre- to star-forming regions, and the 
mid-IR can help to separate these. We find that in the ^=30° field, 
about 55% of the pixels are mIRd, 20% are mIRn, and 25% are 
mIRb, while in the ^=59° field, about 20% are mIRd, 40% are 
mIRn, and 40% are mIRb. There exist significant trends in tem- 
perature and column density between the populations of mIRd to 
mIRb. We find that mIRd dark pixels are about 10 K colder and a 
factor of two higher column density than mIRb pixels. The mIRd 
pixels are likely cold pre- star-forming regions, while the mIRb 
pixels are in regions that have probably begun to form stars. This 
study has shown that Hi-GAL-identified sources span the range 
from cold, pre- star-forming to actively star-forming regions. 

We also include the presence or absence of EGOs (Extended 
Green Objects, indicative of shocks in outflows), CH3OH 
masers, and emission at 8 and 24 jim as star formation tracers. 
We find that warmer pixels are associated with more star forma- 
tion tracers. We also find a wide but plausible trend in tempera- 
ture, where the coldest pixels, on average, are mid-IR-dark, fol- 
lowed by pixels containing EGOs and/or CH3OH masers, then 
8 and 24 fim bright sources. While the systematic errors are too 
large to suggest an evolutionary sequence, this trend is intrigu- 
ing. 

Finally, we identify five candidate far-side IRDCs. These ob- 
jects are cold (~ 20 K), high column density (> 1 x 10^^ cm"^) 
sources identified with Hi-GAL that have very weak or no ab- 
sorption at 8 yum. We explore one such candidate in more detail 
in Section in This object at roughly (^b) = (30.01, -0.27) has 
a high far-IR column density (N(H2) ~ 3 x lO^^cm"^), is cold 
(~ 20 K) and in a position near an abundantly bright mid-IR 



background, yet shows almost no decrement at 8 jim. In fact, 
the 8 Jim extintion-derived column density is almost an order of 
magnitude lower than expected from the far-IR estimate. These 
candidate far-side IRDCs are the first of their kind of which the 
authors are aware. This type of analysis will likely uncover many 
more such objects. With a complete sample of IRDC-like (cold, 
high column density) clouds, independent of the local mid-IR 
background, one could map the clouds in the earliest stages of 
star-formation over the entire Galaxy. 
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